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ABSTRACT 

We present observations of a very massive galaxy at z = 1.82 which show that its morphology, size, velocity 
dispersion and stellar population properties that are fully consistent with those expected for passively evolving 
progenitors of today's giant ellipticals. These findings are based on a deep optical rest-frame spectrum obtained 
with the Multi-Object InfraRed Camera and Spectrograph (MOIRCS) on the Subaru telescope of a high-z 
passive galaxy candidate (pBzK) from the COSMOS field, for which we accurately measure its redshift of 
z = 1 .8230 and obtain an upper limit on its velocity dispersion er* < 326 km s" 1 . By detailed stellar population 
modeling of both the galaxy broad-band SED and the rest-frame optical spectrum we derive a star-formation- 
weighted age and formation redshift of f s f ~ 1-2 Gyr and Zf orm ~ 2.5^4, and a stellar mass of ~ 3^1 x 
1O H M0. This is in agreement with a virial mass limit of M v h < 7 x 1O H M0, derived from the measured cr* 
value and stellar half-light radius, as well as with the dynamical mass limit based on the Jeans equations. In 
contrast with previously reported super-dense passive galaxies at z ~ 2, the present galaxy at z = 1.82 appears 
to have both size and velocity dispersion similar to early-type galaxies in the local Universe with similar stellar 
mass. This suggests that z ~ 2 massive and passive galaxies may exhibit a wide range of properties, then 
possibly following quite different evolutionary histories from z ~ 2 to z = 0. 
Subject headings: galaxies: evolution — galaxies: formation — galaxies: high-redshift 



1. INTRODUCTION 

Understanding the formation of massive elliptical galaxies 
remains a crucial unsolved issue of galaxy evolution. The re- 
cent discovery and the first redshift measurements, through 
deep ultraviolet (UV) rest-frame spectroscopy, of a substan- 
tial population of passively evolving galaxies at z > 1 .4 (e.g. , 
ICimatti et al]l2004t iMcCarthy et d J 120041 IDaddi et al.ll2005l) 
have shown that quenching of star formation in the most mas- 
sive galaxies was already well under way at z — 2. 

A puzzling property of such objects has been revealed 
soon afterwards with some of them being found to have a 
factor of ~ 2 - 5 smaller effective radii compared to local 
early-type galaxies (ETGs) of the same stellar mass (e.g. , 
Daddi et alj|2005t iTruiillo et al.ll2006t lLonghetti et all 120071: 
Cimatti et al.l l2008t Ivan Dokkum et all 120081) . implying that 
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1 Based on data collected at Subaru Telescope, which is operated by the 
National Astronomical Observatory of Japan (S09A-043). 



they are > 10 times denser than their possible descendants in 
the local Universe. Several alternative mechanisms have been 
proposed to make such compact ETGs gro w in size so to fi 
nally meet the propertie s of local ETGs (e.g ., Khochfar & Sill 
2006: iFan et al.l 120081: iNaab et al.l 12009b La Barbera et ~ 
2009; Feldman n et al.l2010l) . but no general consensus has yet 
emerged. 

On the other hand, ETGs at z > 1.4 with large effec- 
tive radii, comparable to the local ETGs, have a lso been 
found (e.g.,lMancini et alj |20!ot ISaracco etal.H 2009. see also 
Dadd i et alj 120051) . indicating a diversity of structural prop- 
erties in the ETG population at z — 2. Moreover, possible 
effects have also been discu ssed that could bia s size estimates 
towards lower valu es (e.g . . IDaddi et al.ll2005|: iHopkins et al.l 
l200^lMancini et al.ll2010tlPannella et al.ll2009l) . 

An independent way to check these issues is by measur- 
ing stellar velocity dispersions (cr*): if high-z ETGs are really 
super-dense, their cr* shoul d be much higher than t hat of local 
ETGs of the same mass. Cappellari et al. (2009) measured 
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er* from deep UV rest-frame spectroscopy of a sample of 9 
ETGs at 1 .4 < z < 2.0 from the GMASS survey dCimatti et alj 
2008), finding two galaxies with similar stellar density and <r* 
as local ETGs of the same mass. The remaining galaxies have 
higher stellar densities and higher cr+ from their stacked spec- 
trum, but still overlapping with the densest local ETGs. 

In this respect, near-infrared (NIR) spectroscopy offers a 
great advantage for the galaxies at z > 1 .4 as the optical break 
at rest-frame 4000 A and Ca II H+K (the strongest spectral 
features of passively evolving galaxies) are redshifted into 
NIR, and the rest-frame optical continuu m is much br ighter 
than t hat in the rest-frame UV. Thus, Ivan Dokkum et al] 

(2009) measured a velocity dispersion of ct* ~ 500 km s -1 
for a compact ETG at z = 2.186 using a deep NIR spectrum 
from|Kriek et al. (2009). This value of ct* is much higher than 
that of the most massive local ETGs, and would be consistent 
with the small half light radius measured for that galaxy. 

In this letter we present a rest-frame optical spectrum of 
a massive, passively evolving high redshift galaxy candidate 
taken with the Multi-Object InfraRed Camera and Spect ro- 
graph (MOIRCS; Echikawa et aLll2006t ISuzuki et al.ll2008l) at 
Subaru Telescope. A cosmology with 51m = 0.3, f^A = 0.7, 
and Ho = 70 km s" 1 Mpc -1 is assumed. 

2. OBSERVATIONS AND DATA REDUCTIONS 

We have obtained 4.7 h ours of MOIRCS spectroscopy of 
34 Bz K-selected galaxie s (IDaddi et al.1 12004) from the cata- 
log of lMcCracken etal.1 (|2010|) in the COSMOS field. We 
used the zJ500 grism with 0"7 slits, providing a resolution 
of R ~ 500 in the /-band over the range 9500-16000A. The 
primary aim of our observation was to measure redshifts for 
passive BzK galaxies (pBzKs) and locate them accurately in 
the COSMOS large scale structure. We preferentially selected 
the most massive pBzKs (which are also the brightest in the 
NIR) to maximize the chance of determining also physical 
information in addition to redshifts for at least a fraction of 
them, in particular the gala xy (#254025) discuss ed in this pa- 
per and already studied by M ancini et al.1 (120101) . The obser- 
vations were made under partly cloudy conditions and with 
i~ l."2 seeing. A sequence of 600s integrations were made 
in a two position dithering pattern separated by 2". The 
AOV-type star HIP 55627 was observed to obtain (relative) 
flux calibration and to correct for telluric absorption. The 
data were reduced using the MCSMDP pipeline (Yoshikawa 
et al., in preparation), including flat-fielding by dome flat, sky 
subtraction between each exposure pairs, bad pixel and cos- 
mic ray rejection, distortion correction, wavelength calibra- 
tion (based on the OH telluric lines), residual sky subtraction 
and finally co-addition with appropriate offsets and weights. 
The 2-dimensional spectra are flux-calibrated using the stan- 
dard star spectrum, and 1 -dimensional spectra were extracted 
with the IRAF apall task using a 1"9 aperture. The abso- 
lute flux calibration was then obtained by normalizing to the 
/-band total magnitude. 

Whilst the results for the complete sample observed with 
MOIRCS will be presented elsewhere we will concentrate 
here on the pBzK galaxy #254025. This galax y one of the 12 
ultra-m assive high redshift ETG candidates in M ancini et al.l 

(2010) , with a photometric redshift of z P hot = 1-71 and very 
bright NIR magnitud es of / AB = 20.32 and K AB = 19.41. 
Manc ini et al.l (120101) also report that the galaxy is non- 
detected at Spitzer/MTPS 24/Ltm to 80 ^Jy implying a star 
formation rate SFR < SOMoyr" 1 . Using HSTI ACS F814W 



imaging (rest-frame UV) they measure a Sersic index n = 
4.1 and an effective radius of r e = 5.7 kpc, consistent 
with the stellar mass -size relation of local elliptical galaxies 
dManciniet al.ll2010h . 

Figure Q] shows the resulting ID and 2D MOIRCS spec- 
tra of #254025. The 4000 A break is clearly seen, together 
with strong Balmer and metallic absorption lines, namely 
Ca II H +He, Ca II K, HS, H 7 , HC, G-band and CN+H0. 
H/3 falls in the region with low atmospheric transmission 
and with strong OH-lines. No emission lines are observed. 
While [O III]AA4959,5007 falls in a region with low atmo- 
spheric transmission, [O II]A3727 is uncontaminated and its 
non-detection sets a 3a upper limit of SFR of ~ 2.5M p 1 yr~ 1 
(not corrected for extinction), using the lKennicutj d 19981) con- 
version. From the spectrum, the absorption line redshift is 
measured as z = 1 .8230 ± 0.0006. 

3. RESULTS 
3.1. Stellar Populations 

Having determined the spectroscopic redshift, we pro- 
ceeded to fit stellar population templates, separately to the 
broad-band spectral energy distribution (SED) and to the 
MOIRCS spectrum. We allowed for a wide range of possi- 
ble star formation histories (SFHs), including: (1) instanta- 
neous bursts, i.e., simple stellar populations (SSP); (2) con- 
stant SFRs for a duration within 0.01-3.5 Gyr 2 , terminated 
by SF quenching and followed by passive evolution; (3) de- 
layed, exponentially declining SFH described as SFR(f,r) oc 
(t /r 2 ) exp(-f /t) with r within 0.01-2 Gyr; (4) exponentially 
increasing SFH, SFR(f,r) oc exp(f/r) for a duration within 
t q = 0.1-3 Gyr, followed by SF quenching and passive evolu- 
tion. We choose r = 0.72 Gyr, corresponding to a stellar mass 
doubling time of ~ 0.5 Gyr, as suggested for z ~ 2 galax- 
ies by the exi stence of tight ste l lar mass-SFR relation, with 
SFR cx- M* (IDaddi et al]|2007t lRenzinill2009l) . For all the 
SFHs, tem plate ages wer e allowed to range in t = 0.4-3.5 Gyr. 
We use a Chabrien d2003l) IMF. The fits were made with metal - 
licities of 0.5xZ Q , Z Q and 2xZ Q . To reduce the number of 
free parameters we have assumed no dust extinction, appro- 
priate for a passively evolving galaxy, considering the strict 
upper limit on the SFR that is set by the spectrum, an assump- 
tion that is validated by the good fit that is achieved in the blue 
continuum (see FigureQ]). 

The SED fitting was carried out for the broad-band 
Biz data from Subaru/Suprime-Cam (Capak et alj 120071 ; 
Tanig uchi et all l2007h. JHK data from CFHT/WIRCAM 
dMcCracken et alJ 120 ICh and the Spitzer /IRAC 3.6pm, 
4.5pm and 5.8pm data (Sanders et al. 2007). Artificial er- 
rors of 0.05 mag for the BizJHK bands and 0. 1 mag for the 
IRAC bands are added to the observed errors in quadrature 
to account for systematics in zero-point determinations, in the 
photometric measurements, and in the stellar population mod- 
els. The templates for the S ED fitting are gen erated fro m pop - 
ulation synthesis models by Maraston (2005, hereafter M05). 

The spectral resolution of the M05 models is significantly 
lower than that of our MOIRCS spectrum . Hence, for fitting 
the spe ctrum, w e used templates from the iBruzual & C hariot 
(I2003L hereafter IBC03I) spectral synthesis library. Although 
these models might not acc ount properly for TP-AGB stars 
(e.g.. iMaraston et aLI 120061) . this effect is not significant at 

2 The age of the Universe at z = 1-82 is about 3.5 Gyr, given the adopted 
cosmology. 
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4000A rest-frame, the wavelength range probed by our spec- 
trum. The template spectra are Gaussian-broadened to an 
overall velocity dispersion of 350 km s -1 (see 33.21 ) to match 
that the observed spectrum, and having fixed it the stellar pop- 
ulation parameters are derived with the \ 2 over the observed 
wavelength range 9500-16000A. 

The stellar population parameters of the best-fit models 
from each adopted SFH are listed in Table Q] and the best- 
fit templates for the spectrum and SED are shown in Figures 
[Hand |2] respectively. The best fit spectra have star-formation 
(SF) weighted ages of t s f = 1 . 14^qj Gyr and stellar masses of 
M* = (2.76^$) x 1O U M for the broad-band SED, and f sf = 
1.88!gj£ Gyr and stellar mass of M* = (3.99^;^) x 10 n M Q 
for the spectrum. The best fit results are from Z = Z Q for 
the SED and Z = 2 x Z Q for the spectrum. However, for the 
spectrum very similar values are derived using solar metal- 
licity models that result in a slightly higher \ 2 (Table 1). We 
notice that solar, or slightly supersolar metallicities are appro- 
priate for local ell iptical galaxies wit h stellar masses similar to 
galaxy #254025 dThomas et al.ll2005l) . The SF weighted age 
of ~ 1-2 Gyr corresponds to an average formation redshift of 
Zf orm ~ 2.5^1, although the SF could have started much ear- 
lier. In the case of the spectral fitting, all SFHs adopted here 
fit equally well with x 2 — 1-3 and f S f ~ 1-2 Gyr, consistent 
with the detection of strong Balmer absorption lines which are 
most prominent in A-type stars. Moreover, in the case of the 
SED fits (apart for the SSP spectra) the various SFHs do not 
give very different \ 2 values and therefore it is not possible 
to discriminate between them. The same can be said for the 
derived metallicities. If we allow for dust reddening, the best 
fitting f S f, and M /Lu would change by only < 10%, with 
some increase of the formal uncertainties within each class of 
SFH. 

3.2. Velocity Dispersion and Dynamical Modeling 

Our high S/N spectrum (~ 8.7 per 60 km s" 1 spectral 
interval in the continuum) allows us to measure a velocity 
dispersion from the absorption line width <7 b s , which is a 
combination of the galaxy stellar velocity dispersion ov and 
the instrumental resolution eri ns t r - This S/N is comparable 
to spectra of G MASS galaxies with su ccessful individual a* 
determinations (Cappella ri et alJ 120091). T herefore we fol- 
lowed the same approach of Cappellari e t all based on the Pe- 
nalize d Pixel-Fitting method (pPXF; Cappellari & Emselleml 
120041). The MILES stellar library containing 985 stars 
dSanch ez-Blazq uez et all 120061) is adopted here since it pro- 
vides the best uniform and complete set of stars. 

Figure [3] shows the best-fit templates from pPXF, corre- 
sponding to cr bs = 350 ± 30 km s" 1 (random) ± 30 km s -1 
(systematic) for the rest-frame wavelength range of 3500- 
4450 A. The random error (la confidence) is determined as 
half of the interval in <7 b s spanned by 68 out of 100 Monte 
Carlo realizations of the input spectrum. A rough estimate of 
the systematic error is obtained as half of the interval in cr b s 
spanned by all repeated extractions of the kinematics using 
different, but equally acceptable, combinations for the values 
of the degree (from 0^1) of the additive and multiplicative 
polynomials in pPXF. Restricting the fit to the region with the 
Balmer and Ca II H+K lines (3700-4100 A in the rest-frame) 
gives cr bs = 300 ± 50 km s" 1 (total error), consistent with the 
value derived from the full spectral range. 

In order to derive cr* we need to determine accurately the 



instrumental resolution. To do this we have used our com- 
bined MOIRCS spectrum without sky subtraction and simul- 
taneously fitted Gaussian profiles to a series of telluric OH- 
lines at A ~ 11,500 A, i.e., near the strongest absorption fea- 
tures of the galaxy's spectru m. The central wavele ngth for 
each OH-line was taken from Rousselot et alJ (120001) and we 
left (Tinsti, OH-line intensities, and constant baseline as free pa- 
rameters. The fitting procedure reproduces the observed sky 
spectra very well. Over A = 3500-4450 A, the instrumental 
resolution changes from 270 km s -1 to 330 km s" 1 . We adopt 
Cinstr = 300 ± 7 km s" 1 (random) ±30 km s -1 (systematic). 
The derived galaxy stellar velocity dispersion is cr* = 

J CTg bs - of nstr , which gives av = 180 ±59 km s" 1 (random) 

±87 km s" 1 (systematic). The relatively large uncertainties in 
(T b s , and as cr b s is close to ai nstr , we cannot place a lower limit 
to cr*. However, we can derive a la upper limit of er* < 326 
km s" 1 (or er* < 385 km s" 1 at the 2a level), which is con- 
sistent with both determinations. If Balmer lines suffer from 
fill-in from emission lines a* could be somewhat lower. 

From the stellar velocity dispersion the virial mass can be 
calculated as M v ; r = Cr e al/G. We have set C = 5 as empiri- 
cally calibrat ed on local galaxies wi th state-of-the-art dynami- 
cal modeling (Cappellari et al. 2006), with the velocity disper- 
sion being measured within a large aperture ( ~ lr e ) as in our 
case. The effective radius r e was measured by Manci nTet alJ 
d201Ch from the HSTl 'ACS F814W image (~ 2900 A in the 
rest-frame) as 0.68" ± 0.07" or 5.7 ±0.6 kpc at z = 1.82. 
Thus the upper limit of the virial mass is derived as M v ; r < 
7.0 x 10 n M Q . 

We have also constructed a dynamical model based on ax- 
isymmetric Jeans dynami cal models as those use d to model 
the GMASS galaxies by ICappellari et all (120091). adoptin g 
a multi-Gaussian expansion (MGE; Emse llem et al.l 1994). 
This method has the advantage that the derived M/L is virtu- 
ally insensitive to possible underestimation of the size, which 
can be a possibility at high redshifts Considering the bolo- 
metric surface brightness dimming of (1 +z) 4 , a factor (1 +z) 
coming from the source redshifting and the Zf-correction be- 
tween rest-frame 2900 A to rest-frame i/-band, we derived 
a rest-frame £/-band luminosity of Lu = 6.7 x 10 u Lo. The 
second moment of the velocity V^. = V 2 + al was also es- 
timated (assuming f3 z = and axisym metry) by usi ng the 
Jeans anisotropic MGE (JAM) method (Cappellari 2008]). The 
upper limit for the dynamical M /Lu can be calculated by 
(M/L[/)j eans = (ov/Vrms) 2 < L0, which can be converted into 
the upper limit of the dynamical mass from the JAM model 
as M Jeans = Lu x (M/Lu )j ean s < 6.8 x 1O H M . Therefore the 
virial mass and Jeans mass agree well though both of them 
are upper limits. A JAM model constructed fr om a noiseless 
model with the best fitting Sersic parameters of Mancini et al. 
(2010), as opposed to the actual ACS image, gives the same 
(M/L[/)j eans within 1%. This is due to the robustness of the 
central M/L recovered using dynamical models (in contrast 
to virial esti mates) to photometric uncertainties at large radii 
(e.g., §3.2 in lCappellari eTaIll2009l) . 

4. DISCUSSION AND CONCLUSIONS 

Figure[4]compares the properties of the galaxy #254025 and 
other z — 2 galaxies for which the same quantities have bee n 
measured ( ICappellari et all l2009t Ivan Dokkum et all 120091) . 
The figure includes ETGs at z — 0.06, selected from the Sloan 
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Digital Sky Survey (SDSS) on the base of their red u—g color 
and high Sersic index n ~ 4 (Blanton et al. 2005). The dynam- 
ical and stellar masses agree very well for the high-z objects, 
within a factor of < 2. Note that our massive galaxy has phys- 
ical properties in good agreement with those of local ETGs of 
similar stellar mass. Our galaxy provides a second example of 
a very massive passively evolving system for which a stellar 
ve locity dispersion has been measured (the GMASS objects 
of ICappellari et al.L which in Figure [4] also lie on the z = 
scaling relations, but have stellar masses below 10 n M Q ). The 
"normal" size and velocity dispersion of our massive ETG is 
strikingly in contrast with the extreme properties (i.e., a very 
high a* = 5 10!^p km s" 1 and a small r R = 0.78 ± 0.17 kpc) of 
the galaxy studied by Ivan Dokkum et al] (120091) with similar 
stellar mass (2 x 1O U M0). This suggests a substantial diver- 
sity in the physical properties of the massive ETG population 
at z ~ 2 including "immature", albeit virialized, systems — 
which will have to evolve into normal z = massive galaxies 
through some physical processes which decrease their veloc- 
ity dispersion and increase their sizes — as well as "mature" 
ETGs, already on the scaling relationships of z = ETGs. It is 
clear that many more observations of similar galaxies are re- 
quired to establish which kind of ETG is commonest at high 
re dshift: either the compac t/high-oy objects like those found 
by Ivan Dokkum et al.l (120091) . or the apparently normal, low- 
er* objects presented in this paper. Also, nothing prevents our 
particular object to evolve further from its present state which 
mimics that of local ellipticals of the same mass. For exam- 
ple, it may grow further and become a brightest cluster galaxy 
(BCG) or a cD galaxy. For this reason, it would be important 
to estimate the volume number density of similar objects at 
high redshifts, and compare it to that of BCGs and cD galax- 
ies. 

To conclude, both very compact ETGs and ETGs following 
the local stellar mass-size and ste llar mass-oy relation s appear 
to co-exist at z > 1.4 (see also iMancini et al.ll2010h . How- 
ever the number of high-z ETGs with individual measurement 
of the velocity dispersion is still extremely small. Increasing 
their sample is of great importance to understand the evolution 
of these galaxies, and in particular how and when they ac- 
quire their final structural and dynamical configuration. This 
paper demonstrates that with reasonable telescope time sev- 
eral absorption features can be detected in the rest-frame op- 
tical spectrum of the high-z ETGs, from which (at least for 
the most massive ETGs) the velocity dispersion and several 
stellar population properties can be derived. 
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Figure 1. MOIRCS spectrum of #254025 at Zspec = 1.82. Top: 4" 2D spectrum; middle: ID spectrum without smoothing (grey line) and with a 25 A binning 
(black thick line); bottom: relative noise level (solid line) and the sky transmission (shaded area). The red solid line shows the best-fit model (see j|3.1> . Positions 
of major emission and absorption lines are indicated by dot-dashed (blue) and dashed (red) lines, respectively, even when not detected. 
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Figure 2. Observed SED for #254025 (blue circles with error bars) compared to the best fit model from Maraston (2005; orange line and symbols). The best fit 
was found for a SFH with delayed exponential star formation rate that continued for about 1 .4 Gyr, with a star-formation timescale of 0. 1 3 Gyr, and for a solar 
metallicity. The parameters of the best fit templates can be found in Table 1 . 
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Figure 3. Results of the pPXF fit for the stellar velocity dispersion of the galaxy. The panels show the resulting fit for the full spectral range (top), and for a 
wavelength range around Ca II H+K (bottom). The black solid line shows the observed spectrum, the red solid line shows the best-fit template, and the green 
diamonds are the residuals (arbitrarily shifted). The blue crosses indicate bad pixels rejected from the fitting. The solid blue line indicates the estimated lcr noise 
level. 
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Figure 4. A comparison of the properties of high-redshift ETGs for which the velocity dispersion has been measured so far (symbols with error bars) with those 
of elliptical galaxies selected from SDSS at z — 0.06 (gray dots), (a): the stellar velocity dispersions vs. the virial masses, (b): comparison between virial and 
stellar masses; the diagonal dashed line corresponds to equality of the two masses, and the dotted lines show a range by a factor of 2. (c): the effective radii vs. 
the virial masses, id): the ef fective radii vs. stellar velo city dispersions. In all panels, the red filled circle represents the galaxy studied here, the green filled circle 
shows the object studied by van Dokkum et al. 1 2009), the blue squares represent the two GMASS galaxies with individual <7* measurements, and the yellow 
triangle represents the properties from the stacked GMASS spectrum (taken from Canpellari et al. 2009). 



Table 1 

Best-fit Stellar Population Parameters for the SED and the Spectrum 
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Note. — Column 1: star-formation history (see 33.lt ; Column 2,8: elapsed time since the onset of star-formation; Column 3,9: star-formation timescale in the case of delayed 
exponential SFH and quenching time for SFHs with constant SFR+quenching and exponentially increasing SFR+quenching; Column 4,10: star-formation weighted age defined by 
f*(T-t)<f>(t) dtj Jq 0(0 dt where T is r on , el and <p{t) is SFR; Column 5,11: stellar mass; Column 6,12: rest-frame f/-band mass-to-light ratio; Column 7,13: reduced \ f° r tne 
best-fit template. 



